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ABSTRACT 


Testing  complete  independence  is  one  of  the  simplest  problems 
onceming  the  covariance  structure  of  a  set  of  measurements.  A 
stepwise  procedure  proposed  by  S.  N.  Roy  and  R.  E.  Bargmann  (1958)  and 
a  trace  criterion  due  to  H.  Nagao  (1973)  are  two  well  known  copetitors 
of  the  likelihood  ratio  test  of  the  hypothesis  derived  assuming  the 
multivariate  normality.  We  consider  some  modifications  of  the  Roy- 
Bargmann  procedure  based  on  combinations  of  independent  tests  and  find 
them  to  be  asymptotically  equivalent  to  the  likelihood  ratio  test,  which 
is  optimal  in  terms  of  the  exact  slopes.  The  operating  characteristics 


of  various  tests  with  samples  of  moderate  size  are  examined  empirically. 
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1.  INTRODUCTION  AND  SUMMARY 

Let  Xj,}^, . ,X^  be  a  random  sample  from  a  p-variate  normal 

population  with  the  covariance  matrix  I.  One  of  the  simplest  problems 
concerning  the  covariance  structure  of  the  multivariate  normal  distribution 
is  of  testing  the  complete  independence  of  the  p  measurements  comprising 
the  vectors  X^.  The  likelihood  ratio  test  for  the  complete  independence 
which  depends  upon  the  determinant  |r|  of  the  correlation  matrix  was 
derived  by  Wilks  (1935) .  The  exact  distribution  of  the  likelohood  ratio 
statistic  is  discussed  and  tabulated  by  Mathai  and  Katiyar  (1979). 

An  alternative  solution  termed  a  step-down  procedure,  which  consists  of 
p-1  independent  tests  was  proposed  by  Roy  and  Bargmann  (1958) .  This 
procedure,  unlike  the  likelihood  ratio  test,  permits  post-hoc  analysis  of 
the  nature  of  dependence  in  case  of  a  rejection  of  the  null  hypothesis 
and  depends  upon  only  the  well  tabulated  F-distribution  for  its  implementation. 

In  this  paper  we  introduce  a  class  of  tests  asymptotically  equivalent, 
in  terms  of  the  exact  Bahadur  slopes,  to  the  likelihood  ratio  test  which  is 
optimal  in  this  sense.  The  presently  available  methods  of  testing 
complete  independence  are  summarized  in  section  2.  The  new  tests  are 
introduced  and  shown  to  be  Bahadur-optimal  in  section  3.  Section  4 
contains  a  Monte  Carlo  comparison  of  these  tests  with  the  likelihood 
ratio  test  and  the  step-down  procedure  when  the  samples  are  of  moderate 
size.  The  empirical  study  also  includes  a  test  proposed  by  Nagao  (1973). 

2.  SOME  TESTS  OF  COMPLETE  INDEPENDENCE 

Let  R  be  the  correlation  matrix  of  a  sample  of  size  N  from  the 
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Np(y,E)  population.  The  likelihood  ratio  test  for  Ho:  £  *  diag  (o2,  a2,  .... 

2 

ffp)  a  rd  rejects  it  if 

A*  |R|'N/2  <  c  ,  (2.1) 

where  the  critical  constant  c  may  be  obtained  from  Mathai  and  Katiyar  (1979) 

or  obtained  by  using  approximations  such  as  Box’s,  and  Bartlett’s  discussed 

by  Mudholkar,  Trivedi,  and  Lin  (1980).  Nagao  (1973)  noted  that  asymptotically 

-2  log  is  a  x  "variable  when  the  null  hypothesis  is  true  but  it  is  an 

asymptotic  normal  variable  for  any  fixed  alternative.  He  suggested 
2  -12 

regarding  x  »  tr(£  Z ^  -I)  which  is  proportional  to  the  variance  of  this 
normal  distribution  as  a  noncentrality  parameter,  i.e.,  a  measure  of 
departure  from  the  null  hypothesis,  and  proposed  a  consistent  estimator 

T  -  ^  tr  (S  SjJ1-!)2  »  ^  tr  (R2-!)  (2.2) 

of  a  multiple  of  t2  as  a  test  statistics  for  H  .  He  obtained  the 
asymptotic  expansion  for  T  in  the  form 

Pr(T  5  x)  ■  pf  +  gTf  ta6  Pf+6  +  a4  Pf+4  * 

a2  Pf+2  +  ao  Pf^  +  °(n"2)  »  (2.3) 

where  f*  p(p-l)/2  and  Pf  »  P[x2f  <  x]  and  afi  *  (p3  -  3p2  +  2p)/12, 

3  2  3 

a^  *  (-2p  ♦  3p  -  p)/8,  and  9.^  ~  (P  ’  P }/4.  He  showed  that  it  is 

satisfactory  for  n  *  100 . 

Roy  and  Bargmann  (1958)  consider  the  null  hypothesis  of  complete 

independence  in  the  form  H„  =  /L  (H  • :  p?  ,,  ,.  ,,  *  0),  where 

o  x«2  01  1  •  1Z....(1-1J 

2 

p^  .  (i-1)  *s  multiple  correlation  between  and  (Xj^, . X^  p  , 

and  note  that  the  sample  step-down  multiple  correlation  coefficients 
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R^  _  j2  ^  are  independently  distributed  when  E  is  diagonal. 

They  propose  rejecting  the  null  hypothesis  when  at  least  one  of  the  component 
hypotheses  is  rejected  by  the  usual  test  for  it,i.e.  when  R?  #  ^  (j.i 
constant.  This  procedure  is  simple  to  implement  as  independently  distributed 


have  variance  ratio  distributions  with  (i-1,  N-i)  d.f.,  i  *  2,  3,  . 


(2.4) 

P  • 


However,  the  procedure  does  require  an  a  priori  ordering  among  the 

measured  variables  and  a  decision  regarding  the  levels  of  the  component 

tests  which,  because  of  the  independence,  are  related  to  allowable  overall 

type  I  error  a  by  (1-a)  *  II  (1-a.).  It  is  common  to  take  a.  «  l-(l-a)^^>”^ 

i=2 

i  »  2,...,p.  Roy  and  Bargmann  gave  the  confidence  bounds  associated  with 

i 

this  step-down  procedure  which  can  be  used  to  gain  an  understanding  of 
the  nature  of  dependence  in  case  Hq  is  rejected. 


3.  A  CLASS  OF  B-OPTIMAL  TESTS 

The  problem  of  the  allocation  of  the  overall  type  I  error  among  the 
component  tests  of  the  step-down  procedure  may  be  avoided  by  considering, 
instead  of  the  variance  ratio  statistics  F^+^,the  P-values  associated 
with  the  individual  tests  i  =  l,...,k,  where  k  *  p-1.  Since  the  statistics 
Fi+1  are  independent  under  H0,  the  P-values  have  independent  uniform 
null  distributions.  These  can  therefore  be  combined  variously  to  construct 
an  overall  test  for  Hq.  The  problem  of  combining  independent  tests  of 
significance  is  classical  and  the  literature  on  the  subject  is  extensive. 

It  is  well  reviewed  in  Liptak  (1958),  Oosterhoff  (1969),  George  (1977)  and 
Mudholkar  and  George  (1979) . 


A  combination  procedure  for  the  P-values  P^,  Pj,...,  P^  associated 

with  k  independent  tests  of  significance  for  hypotheses  H  0^  a  0Q^  vs. 

Hji1  0^  >  0Q^,  i  =  l,2,...,k  is  based  upon  a  combination  statistic 

TCP,,...,?.)  which  is  used  for  testing  the  overall  hypothesis  H  *» (  1 H  .  vs. 

the  alternative  “LJHj^.  The  overall  null  hypothesis  Ho  is  rejected 

when  ¥(P^, . . . .P^)  is  large.  The  following  are  some  of  the  well  known 

combination  statistics:  (i)  The  earliest  proposed  *  min  {-2  log  P^ 

due  to  Tippett,  (ii)  ¥p  =  1-2  log  Pi  due  to  Fisher,  (iii)  ?N  «  *"1(1-Pi), 

*  being  the  c.d.f.  of  standard  normal,  considered  by  Lipkak  (1958), 

and  (iv)  4'l  =  £  log  [P^/Cl-P^)]  introduced  by  George  (1977).  These  statistics 

have  simple  null  distributions.  ¥ ^  is  distributed  as  the  smallest  order 

statistic  of  a  sample  from  the  exponential  population,  ¥p  is  a 

x2  -variable,  has  N(0,  k)  distribution,  and  T,  a  k-fold  convolution  of 
2k  "  L 

logistic  distribution  is  approximately  a  t  variable  with  5k  ♦  4  degrees 
of  freedom.  It  is  easily  seen  that  the  stepdown  procedure  with  equal 
a^'s  is  equivalent  to  the  Tippett  combination  of  its  P-values.  An  account 
of  various  studies  of  the  operating  characteristics  of  combination  methods 
in  the  Neyman-Pearson  and  decision  theoretic  framework  may  be  found  in 
Oosterhoff  (1969) .  However,  none  of  the  methods  can  be  preferred  on  the 
basis  of  these  works.  Littel  and  Folks  (1971)  examined  Bahadur  ARE's 
of  various  methods  and  found  that  among  all  monotone  combination  procedures 
Fisher’s  is  optimal  according  to  this  criterian.  Mudholkar  and  George 
(1979)  showed  that  has  the  same  exact  slope  as  'Fp  and  is  consequently 
optimal.  For  a  recent  account  of  this  aspect  see  Berk  and  Cohen  (1979). 

These  studies  of  asymptotic  relative  efficiencies  concern  combinations 
of  independent  tests;  but  can  be  extended  to  the  methods  of  combining 
tests  which  are  independent  under  the  null  hypothesis  only. 
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The  exact  slope  used  in  defining  Bahadur  ARE  of  a  test  at  an 
alternative  is  the  rate  at  which  -2  log  (P-value  of  the  test)  increases 
with  respect  to  the  sample  size  n,  when  the  alternative  is  true.  Specifically, 
let  large  values  of  a  statistic  Tn  be  significant  in  testing  9  e  |{&  vs. 

Hl:  8  ^  Hq  ,  F0  n(t)  denote  the  distribution  function  of  Tn  and  PR(Tn)  * 

1  -  Fq  n(Tn)  be  the  associated  P-value.  Then  c(0)  *  lim  -2n~*  log  pn(Tn)  , 
when  it  exists,  is  the  exact  slope  of  T  .  c(9)  is  often  obtained  using 

the  following  result  due  Bahadur  (1971,  p.  27). 

Proposition.  Suppose  that  lim  Tn//n  *  b(0)  a.s.  for  each  0  e  H  -  Hq. 

Let  p(t)  *  -n-1log[l-F  (Vn  t)]  and  sappose  that  lim  p  (t)  *  p(t)  exists 

°»n  n-H»  n 

and  is  continuous  on  an  open  interval  containing  the  range  of  b(0).  Then 

the  exact  slope  of  Tr  is  c(0)  »  2  •  p(b(0)). 

Remark.  p(t)  is  sometimes  referred  as  the  index  of  the  sequence 

{T  }  or  of  the  sequence  of  distributions  {F,.  )  of  {T 

n  o,n  n 

Now  consider  the  present  problem  of  testing  the  null  hypothesis  HQ 

that  the  covariance  matrix  £  is  diagonal.  The  step-down  procedure  which 

•  •  2 
involves  testing  the  component  hypothesis  H  . :  p.  .  .  *  0  with 

Ol  1^1 ' lfc • • • 1 

Fi.l  ■  tN-i-lJ/l]  .  [*Vl.l2...1/<1-*2U1.12...i»>  <3» 

i  *  l,2,...,k,  may  be  modified  by  combining  these  independent  (only  under 
H  )  tests  using  a  combination  statistic  ¥(P.  P.  ) ,  where  P.'s  are  the 

W  •  •  •  I  A  X 

P-values  associated  with  F^'s.  We  are  interested  in  the  statistics  of 
the  form  f(P.  P.)  »  £  $  (P.),  where  $*(P.)  *  (P*  )  are  monotone 

*#•••**  i  i  i  i  i)Ti 

decreasing  with  index  p(t)  ■  t,  i  -  1,  2,  ..., 


... 


k.  Let  L  denote  the  family 
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o£  these  tests.  Note  that  if  $  ■  Gjil-P.)  then  under  HQ,  <|i  (P^...,?^) 

is  distributed  as  the  convolution  of  Gj,  G2  . . . ,  G^.  In  fact  L  includes 

in  this  manner  the  combinations  of  the  step-down  tests  based  on  Fisher's 

method  with  G'*(l-t)  *  -2  log  t  and  on  Logit  method  with  G’*(t)  -  log[t/(l-t) ] . 

The  tests  in  the  family  L  are  asymptotically  optimal  and  equivalent 

to  the  likelihood  ratio  test  for  H  .  In  order  to  demonstrate  this,  i.e. 

0 

to  obtain  the  exact  slope  of  iKPj. • • . .P^)  ■  £$(P^),  we  examine  $^(P.) 
which  has  the  same  slope  as  F^+^  given  in  (3.1).  Since  F^+^  is  a  variance- 

ratio  with  (i,  N-i-1)  d.f.,  for  an  alternative  P2^+i.j2  j»  {P^/W-i-l)}1^2 

2  2  1/2 
converges  in  probability  to  f[Pi+1.12  j/Cl-P  i+1>12  *)]/*> 

Moreover,  it  can  be  shown  (eg.,  see  Bahadur  (1971,  p.  13))  that 

-n”1  log(l-F  (»/-nt))  log  (1  +  i  t  ) .  Hence  by  the  above  proposition, 

0  |Tl 

the  exact  slope  of  F^  or  its  monotone  function  4>^(P^)  is 

ci(p2i*l-12...i)  *  'log(1"p2i+1.12...i)  *  C3'2) 

Now,  in  view  of  the  results  by  Berk  and  Cohen  (1979)  it  follows  that  the 

index  of  E#^(Pj)  is  the  same  as  the  index  p^(t)  »  t  of  each  $^(P^)  and 

consequently  the  exact  slope  of  <KP^ . • . • » P^)  is 

cCS-J  -  l»«(l-P21.1.12...i) 

-  -log  1^1  ,  (3.3) 

where  |(P|  denotes  the  determinant  of  the  population  correlation  matrix 


(3.3) 


(P*  (o  — ).  By  particularizing  the  result  in  section  3.4  of  Hsieh  (1979) 
it  is  seen  that  the  likelihood  ratio  test  for  complete  independence  is 
asymptotically  optimal  with  the  exact  slope  -log  j  |,  the  same  as  (3.3). 


4.  AN  EMPIRICAL  EVALUATION 

In  this  section  we  present  a  Monte  Carlo  study  of  the  operating 
characteristics  of  some  of  the  asymptotic  Bahadur  equivalents,  of  the 
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likelihood  ratio  test  for  complete  independence,  described  in  section  3, 

when  the  samples  are  of  moderate  size.  The  study  also  includes 

Nagao's  test  given  in  section  2.  The  finite  sample  behavior  of  the  tests 

are  investigated  in  terms  of  the  power  function  as  well  as  in  terms  of 

the  means  and  s.d.'s  of  the  P-values  of  the  tests  at  various  alternatives. 

The  Monte  Carlo  Experiment.  The  simulation  study  was  conducted 

on  IBM  3032  at  the  University  of  Rochester,  generating  the  random  samples 

from  IMSL  routine  GGNRM.  3000  samples  of  size  n=20  and  n=30  were  drawn 

2  2  2 

from  Np  (Q,  £)  with  *  a2  *****  ap  *  *  and  various  configurations  of 

correlations  from  values  0,  .2,  .4,  .6  and  .8  for  p  *  3,  4,  S.  For  each  sample 

drawn, the  following  test  statistics  were  obtained  from  the  sample  correlation 

matrix  R  =  (r. .) . 

-  ij 

(i)  Likelihood  ratio  based  statistic  t  *  1  -  {N-l-(2p+5)/6}  log  |r|  , 

(ii)  Nagao's  test  statistic  T  *  (N-l)  ^  r?^  , 

(iii)  Step-down  statistics  -  [N-i)/(i-l) ] .  [R  . .  (i-l)^1’^.^  .  .iJ-1^  » 

for  i  *  2,  3,  . .p 

p 

(iv)  Combination  statistic  based  on  Logit  method  ij/L  =  -  ^  logfP^  (1-P^  , 

p 

(v)  Combination  statistic  based  on  Fisher's  method  *  -2  £_  loe  P. 

F  1*2  i  * 

where  P^  are  the  P-values  corresponding  to  the  step  down 
statistic. 

These  test  statistics  were  compared  with  their  critical  constants  determined 
using  the  following  facts  regarding  the  null  distributions: 

(i)  L  is  well  approximated  with  distribution  having  d.f.  p(p-l)/2; 

(ii)  the  critical  constant  Tfl  for  T  may  be  approximated  by 
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T«  -  “  *  S  (  *  <f*2>'f*4» 


2a.  u 

W^2) 


(u+f+2)  + 


2a2U 


2 

where  u  is  the  upper  100a  percentage  point  of  the  x  distribution 

with  d.f.  f  *  p£p-l)/2  ,  and  a2,  a4,  afi  as  given  in  (2.3)  ; 

(iii)  F^  is  distributed  as  a  variance  ratio  F  with  (i-l,N-i) 
d.f. ,  for  i  -  2,  3; 

(iv)  is  approximated  with  a*t^  ,  a  constant  times  student's  t, 

where  a  =  v  {k(5k*2)/(3(Sk+4)}1/2  ,  k  *  p-1,  and  d.f.  v  =  5k+4  ; 

2 

(v)  <|/p  is  distributed  as  a  x  with  2 (p-1)  d.f. 

The  pow  of  each  test  was  estimated  by  the  proportion  of  times  the 

null  hypothesis  was  rejected  by  the  corresponding  test.  The  s.d.  of  any  of 

-1/2 

these  estimates  I  {(3000)x4>  .  The  P-values  corresponding  to  the 

tests  were  obtained  using  the  equation  (2.3)  and  the  results  on  the  null 
distributions  as  mentioned  above.  The  P-values  in  each  case  were  averaged 
and  their  standard  deviation  was  computed. 

The  estimated  power  functions  and  the  means  of  the  P-values  of  the 
five  tests  at  various  alternatives  are  given  in  Table  1  and  Table  2 
respectively.  The  Monte  Carlo  experiment  with  3000  simulation  was  first 
conducted  with  n=20  and  p*3  for  the  correlation  configurations  appearing 
in  the  tables.  After  an  examination  of  the  results  it  was  performed  with 
p*4  and  5  for  two  special  configurations,  namely  (i)  the  extreme  configuration 
in  which  only  the  first  correlation  configuration  p2^  is  nonzero,  and 
(ii)  the  symmetrical  configuration  where  all  correlation  coefficients  are 
equal.  As  a  confirmation  of  the  findings,  the  procedure  was  repeated 


with  n=30.  The  s.d.'s  of  the  P-values  and  the  results  for  n=30  which  are 
not  included  in  this  paper  are  available  from  the  authors. 

Conclusions.  Two  features  of  the  comparative  behavior  of  the  five 
tests  clearly  emerge  from  the  two  tables:  (i)  In  case  of  the  extreme 
configuration,  with  pjj  3*  0,  p^  *  0  otherwise,  the  step  down  procedure 
is  preferable.  Its  superiority  over  the  other  four  tests  increases  as 
p  increases.  Nagao's  test  is  the  poorest  in  this  case,  (ii)  Nagao's 
test  dominates  others  if  HQ  is  violated  in  a  symmetric  manner,  i.e.,  when 
p„'s  are  nonzero  and  equal.  The  stepdown  test  is  the  weakest  in  this  case. 
It  is  also  observed  that  the  likelihood  ratio  test,  and  the  two  combinations 
of  the  P-values  of  the  stepdown  components  are  generally  comparable  and 
are  preferable  except  against  the  two  special  alternatives. 
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Table  1.  The  Empirical  Power  Functions  for 
Samples  of  Size  n=20  with  Monte  Carlo  of  Size  3000 


Nonzero 

L.R. 

Nagao 

Step  Dn 

Logit 

Fisher 

Corr. 

Test 

Test 

Test 

Comb. 

Comb. 

Ho 

.0490 

.0480 

.0527 

.0503 

.0473 

»21  " 

■  .4 

.088 0 

.0900 

.0993 

.1003 

.1027 

.2673 

.2650 

.3300 

.2957 

.3227 

*  .6 

.6503 

.6357 

.7540 

.6437 

.7217 

»  .8 

.9823 

.9780 

.9930 

.9717 

.9877 

P31  * 

*  .4 

.0963 

.0957 

.0867 

.0893 

.0883 

.2710 

.2663 

.2420 

.2197 

.2347 

*  .6 

.6733 

.6493 

.6470 

.5347 

.6107 

*  .8 

.9800 

.9770 

.9757 

.9390 

.9717 

p32  m  *? 

»  .4 

.0990 

.0967 

.0930 

.0877 

.0890 

.2713 

.2660 

.2373 

.2280 

.2350 

*  .6 

.6570 

.6430 

.6273 

.5393 

.5953 

*  .8 

.9717 

.9650 

.9707 

.9210 

.9587 

C21  '  °31  ■ 

•  .4 

.1393 

.1417 

.1320 

.1420 

.1403 

.5983 

.5400 

.5270 

.6200 

.6170 

*  .6 

.9990 

.9930 

.9933 

.9977 

.9990 

“si  ■  "32  ■  -J 
■  .4 

.1387 

.1357 

.1187 

.1173 

.1220 

.5967 

.5403 

.5663 

.4767 

.5407 

*  .6 

.9980 

.994  7 

.9977 

.9817 

.9973 

p2i  -  p32  * 

*  .4 

.1583 

.1550 

.1543 

.1660 

.1673 

.5843 

.5410 

.5203 

.6050 

.5987 

■  .6 

.9987 

.9940 

.9923 

.9980 

.9987 

All  p . . ’ s  =  .2 

.1890 

.2150 

.1693 

.1933 

.1937 

^  *  .4 

.6340 

.6847 

.5353 

.6453 

.6353 

*  .6 

.9490 

.9650 

.9073 

.9567 

.9527 

*  .8 

1.000 

.9997 

.9993 

.9997 

.9997 

H 

0 

.0473 

.0457 

.0457 

.0480 

.0463 

°21  ' 

*  .4 

.0763 

.0803 

.0960 

.0953 

.1020 

.2013 

.1937 

.2817 

.2350 

.2707 

*  .6 

.5227 

.4837 

.7150 

.5530 

.6643 

■  .8 

.9400 

1 _ 

.8993 

.9840 

.9167 

.9727 

Continued 
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Table  1  Continued 


p 

Nonzero 

Corr. 

L.R. 

Test 

Nagao 

Test 

Step  Dn 
Test 

Logit 

Comb. 

Fisher 

Comb. 

4 

All  p..'s  = 

.2 

.2253 

.2860 

.1680 

.2283 

.2257 

.4 

.7270 

.8153 

.5823 

.7440 

.7343 

M 

.6 

.9833 

.9923 

.9427 

.9860 

.9837 

* 

.8 

1.000 

1.000 

.9997 

1.000 

1.000 

5 

H 

o 

.0470 

.0463 

.0447 

.0527 

.0473 

P21  * 
a  .4 

.0670 

.0677 

.0897 

.0817 

.0857 

.1520 

.1493 

.2517 

.2043 

.2313 

a  .6 

.3923 

.3607 

.6760 

.4573 

.S823 

a  .8 

.8727 

.7377 

.9840 

.8673 

.9597 

All  p . . ' S  a 

.2 

.2773 

.3820 

.1887 

.2867 

.2763 

1J  • 

.4 

.8223 

.9080 

.6150 

.8347 

.8220 

as 

.6 

.9897 

.9963 

.9537 

.9903 

.9897 

as 

.8 

1.000 

1.000 

1.000 

1.000 

1.000 

Table  2.  Estimated  Means  of  the  P-values  at  Various  Alternatives  for 
Samples  of  Size  20  with  Monte  Catlo  of  Size  3000 


Nonzero 

L.R. 

Nagao 

Step  Dn 

Logit 

Fisher 

p 

Corr. 

Test 

Test 

Test 

Comb. 

Comb. 

3 

Ho 

.4969 

.4975 

.4978 

.4956 

.4971 

p» :  :l 

.4345 

.4340 

.4250 

.4278 

.4239 

.2482 

.2485 

.2201 

.2418 

.2213 

*  .6 

.0775 

.0805 

.0555 

.0863 

.0610 

*  .8 

.0048 

.0062 

.0023 

.0069 

.0033 

P31  *  '? 

■  .4 

.4229 

.4215 

.4314 

.4339 

.4302 

.2544 

.2553 

.2750 

.2926 

.2749 

«  .6 

.0793 

.0823 

.0900 

.1234 

.0968 

*  .8 

.0048 

.0061 

.0052 

.0155 

.0070 

p32  * 

*  .4 

.4221 

.4216 

.4334 

.4330 

.4310 

.2531 

.2543 

.2739 

.2916 

.2736 

■  .6 

.0834 

.0832 

.0919 

.1248 

.0975 

-  .8 

.0059 

.0076 

.0065 

.0170 

.0087 

“21  *  P31  - 

*  .4 

.3659 

.3658 

.3656 

.3708 

.3610 

.0984 

.1081 

.1117 

.1000 

.0950 

«  .6 

.0007 

.0035 

.0028 

.0008 

.0007 

P31  ’  p32  ' 

•  .4 

.3627 

.3641 

.3791 

.3864 

.3773 

.1009 

.1102 

.1141 

.1459 

.1190 

■  .6 

.0007 

.0034 

.0009 

.0069 

.0012 

P21  =  P32  =  *? 

*  .4 

.3606 

.3613 

.3607 

.3637 

.3562 

.1011 

.1101 

.1134 

.1018 

.0976 

*  .6 

.0007 

.0034 

.0027 

.0008 

.0006 

All  p .  . 1  s  =  .2 

.3217 

.3152 

.3340 

.3250 

.3237 

*  .4 

.0985 

.0871 

.1183 

.0970 

.0978 

*  .6 

.0115 

.0084 

.0192 

.0112 

.0114 

«  .8 

.0001 

.0000 

.0003 

.0001 

.0001 

4 

Ho 

.5030 

.5017 

.5037 

.5064 

.5047 

P21  »  .2 
■  .4 

.4509 

.4491 

-.4291 

.4405 

.4312 

.3090 

.3105 

.2523 

.2929 

.2586 

*  .6 

.1202 

.1280 

.0689 

.1191 

.0812 

*  .8 

.0124 

.0204 

.0031 

.0177 

.0057 

Continued 
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Table  2  Continued 


B 

Nonzero 

Con. 

L.R. 

Test 

Nagao 

Test 

Step  On 
Test 

Logit 

Comb. 

Fisher 

Comb. 

H 

All  p.  »s  =  .2 

.2876 

.2705 

.3145 

.2910 

.2904 

-  .4 

.0660 

.0484 

.0981 

.0643 

.0653 

■  .6 

.0040 

.0018 

.0125 

.0034 

.0039 

*  .8 

.0000 

.0000 

.0002 

.0000 

.0000 

5 

Ho 

.5025 

.5047 

.5012 

.5035 

.5029 

P21  *  *2 

.4640 

.4643 

.4400 

.4483 

.4397 

*  .4 

.3355 

.3346 

.2626 

.3056 

.2724 

«  .6 

.1658 

.1769 

.0838 

.1534 

.1008 

«  .8 

.0251 

.0451 

.0041 

.0288 

.0089 

All  p..'s  »  .2 

.2639 

.2334 

.3042 

.2674 

.2668 

^  *  .4 

.0407 

.0236 

.0832 

.0393 

.0414 

*  .6 

.0023 

.0008 

.0099 

.0018 

.0020 

*  .8 

.0000 

.0000 

.0002 

.0000 

.0000 
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